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In this paper we consider the effects of adding curvature in extended cosmologies involving a 
free-to-vary neutrino sector and different parametrizations of Dark Energy (DE). We make use of 
the Planck 2018 cosmic microwave background temperature and polarization data, Baryon Acoustic 
Oscillations and Pantheon type Ia Supernovae data. Our main result is that a non-flat Universe 
cannot be discarded in light of the current astronomical data, because we find an indication for a 
closed Universe in most of the DE cosmologies explored in this work. On the other hand, forcing 
the Universe to be flat can significantly bias the constraints on the equation of state of the DE 
component and its dynamical nature. 
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I. INTRODUCTION 


On large scales, our Universe is almost homogeneous 
and isotropic, but as far as its curvature is concerned we 
cannot firmly conclude that it is spatially flat [1]. Even 
though the observational probes in the past years were 
in agreement with spatial flatness [2-7], however, this re- 
sult has recently been questioned by some of the exper- 
iments. The investigations with the Cosmic Microwave 
Background (CMB) temperature and polarization spec- 
tra from Planck 2018 team using the baseline Plik likeli- 
hood are suggesting that our Universe could have a closed 


geometry at more than three standard deviations [8- 
11] (40 when the physical curvature density is consid- 
ered [12]), and such indications are coming mainly from 
the temperature data affected by an otherwise inexplica- 
ble excess of lensing (Alens problem [11, 13, 14]). Further- 
more, also from the complementary examinations using 
the alternative CamSpec likelihood [15, 16], the closed 
geometry of the Universe is supported at more than 99% 
CL, and the indication seems to persist even consider- 
ing the CMB temperature-only data of the new Planck 
PR4 analysis [17]. Additionally, an indication for a closed 
universe is also present in the BAO data, using Effective 
Field Theories of Large Scale Structure, once the (hid- 
den) assumptions of flatness (in the fiducial cosmology, in 
the reconstruction process, and in the covariance matrix) 
are removed from the beginning [18]. These outcomes 
challenged the assumptions of flatness present in the 
standard ACDM model and sparked a debate about the 
flatness of our Universe. In fact, if there is no solid the- 
oretical argument, then observations in agreement with 
a small curvature cannot be used as a proof for a spa- 
tially flat universe [19], and the assumption of a spatially 
flat Universe on the spatially curved Universe may sig- 
nificantly affect the cosmological parameters [20]. There- 
fore, in this article we want to investigate, in a systematic 
way, some well known and classical cosmological models, 
leaving the curvature of the Universe as a free parame- 
ter and allow the observational data to pick up the best 
possibilities, removing possible biases due to the flatness 
assumption. Moreover, we are also interested in inves- 
tigating what happens to the current cosmological ten- 
sions [21] when the curvature of the Universe is allowed 
into the picture (also see [22-30]). In particular, we want 
to see what happens to the value of the Hubble constant 
and the > 5c tension existing between the ACDM based 
Planck 2018 estimate [8] and the SHOES (Supernovae and 
Ho for the Equation of State of dark energy) measure- 
ment [31, 32], known as the Hubble tension [33, 34] (see 
also Refs. [35-69] aiming to alleviate the Ho tension in 
various ways), together with the > 30 tension in the Ss 
parameter [70] defined as à combination of the ampli- 
tude of the matter power spectrum og with the matter 
density at present Qmo (Ss = a8 y Nmo/0.3) between the 
CMB data (within the ACDM assumption) and the weak 
lensing experiments [71-74] (see Refs. [40, 54, 63, 75-83] 
offering a possible route to alleviate the Sg tension in 
various alternatives to ACDM). 


Following this approach, in this article, we have consid- 
ered a variety of cosmological models in the curved back- 
ground (i.e. allowing the curvature of the Universe as a 
free parameter), namely, the standard ACDM, wCDM, 
the Chevallier-Polarski-Linder parametrization [84, 85], 
where the dark energy has a dynamical equation of state 
parameter, and lastly a recently introduced emergent 
dark energy model, known as the phenomenologically 


emergent dark energy model [86, 87] !, and their vari- 
ous extensions by allowing a free-to-vary neutrino sector 
characterized by the total neutrino mass and the effec- 
tive number of relativistic degrees of freedom. The mod- 
els, together with their extensions in the neutrino sector, 
have been constrained consistently using the CMB data 
from Planck 2018, Baryon Acoustic Oscillations distance 
measurements, and the Pantheon sample of the Type Ia 
Supernovae. This work is complementary to similar ap- 
proaches used in Refs. [11, 14, 54, 103], but differs in 
the light of combination of parameters explored, obser- 
vational datasets, and the extended analysis in the dark 
energy sector. In particular, for the first time, we are in- 
cluding in the analysis a dynamical dark energy equation 
of state together with a non-zero curvature parameter 
that is free to vary. 

The paper has been structured as follows. In sec- 
tion II, we introduce the framework we are exploring; 
in section III we describe the observational data and the 
method used for the statistical analysis; in section IV we 
present our main results. Eventually, we end in section V 
with a summary of them and our Conclusions. 


II. DARK ENERGY IN A CURVED UNIVERSE 


As argued in the "introduction", we start with the ho- 
mogeneous and isotropic spacetime characterized by the 
Friedmann-Lemaitre-Robertson-Walker (FLRW) line el- 
ement given by 


dr? 


ds? = —dt? + a?(t) |  —— 
s +a (t) EEE 


+ r?(d6? + sin? 6d¢*)| , 
(1) 
where a(t) is the expansion scale factor of the uni- 
verse and K denotes its curvature scalar. The values 
K = 0, 41, —1, respectively correspond to a spatially flat, 
closed and an open geometry of the universe. In the con- 
text of General Relativity, we assume that the matter 
distribution of the universe is minimally coupled to the 
gravity and none of the fluids are interacting with any 
component of it. Then the dynamics of the universe can 
be described from the Einstein's gravitational equations 


K 8rG 
H? =-— + i 2 
: K 
2H --3H? = — - i 
+3 "OE 8TG Sop 35 (3) 


where an overhead dot represents the time derivative; 
H = a(t)/a(t) is the Hubble function of the FLRW uni- 
verse; p; and p; are, respectively, the energy density 


1 Apart from these dark energy parametrizations, a large variants 
of works aiming to probe the dynamical nature of the dark energy 
are available in the literature, see for instance Refs. [88-102]. 


and pressure term of the i-th component; G is the New- 
ton’s gravitational constant. Now introducing the criti- 
cal density, pe = 3H?/87G, eqn. (2) can be expressed as 
1 = Qg +}; Qi, where Q; = pi/ pec is termed as the den- 
sity parameter for the i-th fluid and Qg = —K/(a? H?) is 
termed as the curvature density parameter.? Notice that 
the evolution of the energy density of each component 
can be found by solving the balance equation, 


dpi + 3H 2 pit p) = 0, (4) 


which can be obtained from the field equations (2) and 
(3). Since the fluids are not interacting with each other, 
therefore, we have 6; + 3H (p; + pi) = 0 for each i. 
Here, we consider that the total energy density of the 
Universe is distributed among radiation (r), baryons (b), 
neutrinos (v), cold dark matter (CDM) and a dark energy 
(DE) component. For i = r,b,v, CDM, DE, one can find 
the evolution law of the energy density by solving the 
conservation equation for the i-th fluid. Now, concerning 
the dark energy sector, using the conservation equation 
(4) for its equation-of-state wpe = ppg/pprg, one can 
find the evolution of the energy density as 


“1+ wol’) gy 
12 


PDE = PDE,0 €Xp E I ? (5) 
0 

where ppg,o is the present value of the dark energy den- 
sity. Notice that for wpg — —1, the energy density be- 
comes constant which represents a cosmological constant 
in the Universe. As the nature of dark energy is not 
yet clearly understood to anyone, thus, one might be in- 
terested to see how different variants of the dark energy 
model can affect the evolution of the expansion history 
along with the curvature of the Universe. As there are 
many proposals of dark energy in the literature, here, we 
pick up some well known choices as follows: 


e The cosmological constant as the dark energy can- 
didate characterized by its equation of state wpg — 
—]. 


e The dark energy with a constant equation of state 


other than the cosmological constant (i.e. wor Æ 
-1)3 


e The case when the dark energy equation of state is 
dynamical and it assumes the Chevallier-Polarski- 
Linder parametrization [84, 85]: 


WDE = wo + Wa(l — a) , (6) 


2 At this point we mention that unless otherwise specified, “0” 
attached to any quantity will refer to its present value; for ex- 
ample, Q go refers to the present value of the curvature density 
parameter. 

3 Let us comment that from now on we shall use w instead of wpg 
and we label the corresponding cosmological model as wCDM 
(instead of wpgCDM) as already done in the literature. 


where wo is the present value of wpg and wa = 
d'wpg/da at a = 1 is another free parameter in this 
parametrization.^ 


e Finally, we consider a very exotic dark energy 
model, namely the phenomenologically emergent 
dark energy (PEDE) where DE does not have any 
effective presence in the past but it emerges at late 
times. In this model scenario, the dark energy den- 
sity (ppg) is directly parametrized which takes the 


form [86, 87] 
i = QDE 0 [1 — tanh(logio(1 + z))| , (7) 
c,0 


where pec,o is the critical energy density at present 
time and Qpr o is the present value of the DE den- 
sity parameter Npe = ppE/Perit- We note that even 
if this parametrization has a dynamical nature, this 
does not offer any additional degree of freedom, or 
free parameter. 


Having all the models, one can now constrain them in 
presence of the curvature of the Universe. In the next 
section we describe the observational datasets and the 
statistical process aimed to confront the models with the 
data. 


IIl. OBSERVATIONAL DATA AND 
STATISTICAL PROCESS 


In this section we describe the main observational 
datasets that have been used to constrain the underlying 
cosmological models and the methodology of the statis- 
tical simulations. In what follows we first describe the 
observational datasets. 


e Cosmic Microwave Background (CMB) Ob- 
servations: We consider the measurements of the 
CMB observations from the final release of the 
Planck 2018 team. In particular, we have used the 
CMB temperature and polarization power spectra, 
that means plikTTTEEE-+lowl+lowE from Planck 
2018 [8, 104]. 


e Baryon Acoustic Oscillations (BAO): Various 
measurements of the BAO data have been used in 
this work. In particular we have used the measure- 
ments from 6dFGS [105], SDSS-MGS [106], and 
BOSS DR12 [107] as considered by the Planck 2018 
team [8]. 


4 Note that the corresponding cosmological model is labeled as 
wowa CDM. 


e Pantheon sample from Supernovae Type Ia 
(SNIa): Along with the CMB and BAO measure- 
ments, another potential cosmological probe is the 
SNIa, the fist observational data indicating the ac- 
celerating expansion of the Universe. Here we took 
a recent collection of SNIa, known as the Pan- 
theon sample [108] distributed in the redshift in- 
terval z € [0.01, 2.3]. 


Now we discuss the parameter space of each cosmolog- 
ical scenario considered in this article. We have consid- 
ered the extensions of ACDM, wCDM, wow; CDM and 
PEDE by allowing the curvature of the Universe (Qg) 
and a free-to-vary neutrino sector. The free-to-vary neu- 
trino sector has two free parameters: M, = X` m, (the 
total neutrino mass) and Neg (the effective number of 
relativistic degrees of freedom ). The scenarios ACDM 
+ Ox, ACDM + Ox + M,, ACDM + Og + Nog and 
ACDM + Ox + M, + Neg, respectively contain 7, 8, 8, 
and 9 free parameters as follows: 


£2 { Oh? Qh? 1008546, T, ns, log[109 A,], Qx}, 
Ly = L1 U {M,}, 
L3, = Li U Neg}; 
La = Lı U {M,} U {Neg}. 
The scenarios CDM + Qg, wCDM + Qg + M,, 
wCDM + Qg + Neg and wCDM + Or + My + Neg, 


respectively contain 8, 9, 9, and 10 free parameters as 
follows: 


w= { ph, Q2, 1000.46, T, ns, log[101 A,], Qe, wh, 


Wz = W1 U {M,}, 

W3,= W1 U {Nef}, 

W4 = W1 U {M,} U {Neg}. 
The scenarios wowaCDM + Or, wowaCDM + Qg + 
M,, wowa CDM + Or + Neg and wowa CDM + Or ep 


My + Neg, respectively contain 9, 10, 10, and 11 free 
parameters as follows: 


Ci = [ou QU. 10006, T, ns, log[10'^ A,], Ox, Wo, wa, 


C9 =C U (M, 

C3, = CU {Neg}, 

C=C U iM, U {Nes}, 

Finally, the scenarios, PEDE + Og, PEDE + Or + Mr, 
PEDE + Ox + Neg and PEDE + Qx + My, + Neg, 


respectively contain 7, 8, 8, and 9 free parameters as 
follows: 


&= {05h?, Qh, 100050, T, ns, log[10'? A,], Qx}, 
E2 = &4U {Mp}, 

E3, zéQU {Nef}, 

ELZ E&U {M,} U {Neg}. 


To constrain the parameter spaces of the cosmologi- 
cal scenarios described above, we have used the publicly 
available Markov Chain Monte Carlo code CosmoMC [109, 
110] (see http://cosmologist.info/cosmomc/), which 
supports the Planck 2018 likelihood [104] and addition- 
ally has a convergence diagnostic following the Gelman- 
Rubin statistics [111]. We adopt on the parameters dis- 
cussed above the flat priors listed in Table I. 


Parameter Prior 
Oph? [0.005 , 0.1] 
Qh? |[0.001, 0.99] 

100 6c [0.5 , 10] 
T [0.01 , 0.8] 
log(107° As) | [1.61, 3.91] 
ns [0.8 , 1.2] 
Ox [-0.3 , 0.3] 

wo [-3, 0] 

Ua [-3, 3] 

M, [0 , 1] 

Neg [2.2, 4] 


TABLE I. List of the flat priors assumed on the independent 
parameters. 


IV. OBSERVATIONAL CONSTRAINTS 


This section is entirely devoted to describing the ob- 
servational constraints on various dark energy scenarios 
in a curved Universe. In what follows we report the con- 
straints on each model and its various extensions in a 
dedicated way. 


A. Non-flat ACDM and its extensions 


In this section we describe the various extensions of 
non-flat ACDM scenario. We note that for all extensions 
we have made use of the same datasets and their combi- 
nations and for each parameter either we report its 68% 
and 95% CL constraints or the upper/lower limits at 68% 
and 95% CL. 


1. ACDM + Ox 


In Table II we show the constraints on the ACDM + 
Qx scenario for various datasets, and in Figure 1 the 1D 
and 2D contour plots for a few important parameters. 


Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 
OT RNEMEEYTCUDOHID  0.1197+0 0014+029 0119140001540007 NDERIT REESE: 
Qu? —_0,022619.00037+9.00082 g.qragqt9.00017+9.00020 9, qp2qs+0-20018+0-000%2 0,0294g 0-000150 00030 
006  1.04117+2:00034+0:00065 1 ggggs+0-200%2+0-00082 1 gqiq3+0-000%2+0-00081 ] oqgg+?,0081+0. 90082 
: Ju MET 48-0076+0.010 MET I MT ru 
n. ru MEM EI MEC u MELLE 
EE MEOS MEET EEUU. 
2x0 —0.043* 5 15.034 0.000870 oo19.-0.0039  —0-00607 5 0542 0.012 0-008 fy po19.- 0.0038 
Omo 0.4811 0.067 -0.12 0.3096 + 0.0066—0.012 0.33670 024.-0.041 0.3091000680012 
s dn c ESO — E o ES s 
Ho [Km/s/Mpc] 54.5* $97 72 67.90 67 13 65.2155 14 67.97 6521.3 
Ss 0.978 5 45 0.096 0.82475 514. 0.024 0.8503 $025 o. o47 0.8233 5 012 0.024 
147.34 9-31-4-0.60 147.16+0-33+0.58 147.24+0-31+0.61 


Tdrag [Mpc] —0.31—0.61 


—0.30—0.62 


147.177 9:30x0.60 


—0.31—0.61 0.30—0.59 


TABLE II. 68% and 95% CL constraints on various free and derived parameters of the ACDM + Qx scenario for several 


observational datasets. 


We start investigating the constraints from the CMB 
alone case which are shown in the second column of Ta- 
ble II. We see that for CMB alone, a strong indica- 
tion of a closed Universe is found at more than 95% CL 
(Oxo = —0.043*00:1). The Hubble constant takes a 
very low value (Ho = 54.5733 km/s/Mpc at 68% CL) 
increasing Hubble tension, and due to the strong correla- 
tion between Ho, Qg and Qmo, the matter density Qmo 
takes a higher value as expected. Also we notice that the 
tension in Sg increases for this case. So, effectively, we 
see that CMB alone gives a strong indication for a closed 
Universe at the expense of increased tensions in both Ho 
and Sp. 

When BAO data are added to CMB, that means when 
we consider the combination CMB+BAO (see the third 
column of Table II), we see that the Qg becomes consis- 
tent to the spatially flat Universe and the Hubble con- 
stant goes up leading to Planck's ACDM value. However, 
this result is obtained with the combination of datasets 
in disagreement at more than 3c [9, 10, 112], making 
the combination CMB--BAO not safe and the results are 
not completely reliable. The results for CMB+Pantheon 
and CMB+BAO-+Pantheon remain almost similar un- 
like a mild lower value of the Hubble constant in the 
CMB-+ Pantheon case (Ho = 65.277} km/s/Mpc at 68% 
CL for CMB+Pantheon) but in both the cases, the spa- 
tial flatness is found to be consistent with these data. 


2. ACDM+QxK + M, 


The first non-flat ACDM extended scenario that we an- 
alyze involves massive neutrinos. We refer to this model 
as ACDM + Qg + M,. In Table III we provide the ob- 
servational constraints on the parameters, while in Fig- 


ure 2 we show their 1D and 2D marginalized posterior 
distributions. 


Also in this extended model, considering only 
the Planck satellite data on the CMB temperature 
anisotropies and polarization, we obtain a preference for 
a curved cosmological spacetime at more than 9596 CL, 
with the constraint on the curvature parameter reading 
Oxo = —0.073*0099 at 9596 CL. As for the previous 
baseline case, this preference is strongly reduced when 
the Planck data are combined with the BAO datasets 
that are in tension with it, and their combination prefers 
spatial flatness within one standard deviation. As con- 
cerns the value of the expansion rate today, due to its 
degeneracy with the mass of the active neutrinos, in this 
extended parameter space the CMB measurements pre- 
fer further lower values of Ho giving a less tight bound 
Ho = 48.3725 km/s/Mpc at 6896 CL, which is also in 
strong tension with the SHOES independent measure- 
ment of Ho [31], at the level of ~ 4.30. Combining the 
CMB measurements with the BAO and Pantheon data, 
the tension between the different datasets remains statis- 
tically significant. In particular the values inferred com- 
bining CMB and BAO, Hyg 67.84 + 0.67 km/s/Mpc 
at 68% CL, is similar to the result obtained in the flat 
ACDM model, while combining CMB and Pantheon we 
get Ho = 65.3153 km/s/Mpc at 68% CL. Notice that 
both these values are in tension with the SHOES collab- 
oration at the level of 4.30 and 3.20, respectively. Fur- 
thermore, due to the strong anti-correlation between Ho 
and the matter density parameter, lower values of Ho 
prefer higher values of Qmo, resulting in a severe tension 
also for the Sg parameter. Comparing the bound of Sg 
derived by the CMB measurements alone with respect to 
the constraint obtained including also the BAO and Pan- 
theon data, we can observe a tension at the level of about 
^ 3c. Finally, concerning the neutrino masses, the CMB 
data alone constrain their total mass to be M, < 0.79 
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FIG. 1. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario ACDM + Qx for various observational datasets. 


eV at 95% CL, i.e. a factor more 3 relaxed than the 
flat ACDM+M, scenario However, neutrinos suppress 
structure formation at small scales and therefore astro- 
physical galaxy clustering measurements turn out to be 
crucial to improve the constraining power on their to- 
tal mass. For this reason, the upper bounds on M, can 
be significantly improved to M, < 0.17 eV (always at 
95% CL) including also the Baryon Acoustic Oscillation 
and the Pantheon datasets. This upper limit should be 
compared with that obtained for a flat ACDM+M, sce- 
nario, that is M, < 0.13 eV at 95% CL for CMB+BAO. 
In other words, the assumption of flatness produces a 
much stronger upper limit on the total neutrino mass 
than when the curvature can vary, biasing conclusions 
that are important to laboratory experiments. 


The second non-flat extended cosmological model an- 
alyzed here, accounts for the possibility to have a larger 
effective number of relativistic degrees of freedom Neff 
at recombination. We recall that Neg is defined by the 


relation 
7/4\4/3 

Prad = Nets (=) Py (8) 
with p, the present Cosmic Microwave Background en- 
ergy density. Within the standard ACDM cosmological 
model, this parameter is fixed to Neg = 3.044 [113-118] 
which consists of three massless neutrino species and an 
additional correction coming from the non-instantaneous 
neutrino decoupling. Here we regard Neg as an addi- 
tional free parameter of the ACDM + Qg + Nog cos- 
mological model. It should be noted that when testing 


Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 
a a a E — QI00 8 DEDOS Qro 001440, 0089 
Qu  0.02254+0:20018+0.00038 g.qragqt9.00018+9.00020 0.0294gF0.00018+0.00032 9, gpnq1+0,00018+0.000%1 

1008ra  1.04101+0:00034+0:00068 ggg +0:00031+0:00061, j 9413+0:00032+0:00063 ] qggg7+0,00088+0.00088 
ZEND MD a MED EE 
ns 0.969375 0050.-0.0099  0-9659*5.0045.0.0057 0-960787 5.0048-0.008s 0.90637 o Qo46--0. 0088 

EU S S o S 
Qo —0.073* 5 023 0.071 0.0007 0 0o22.-0.0041  —0-0061* 5 0056-0.012 0.0007 5 o024.-0.0045 
Omo 0.65* 531 -0.28 0.3102 0 oo69 0.014 0.334* 0023. 0.042 0.30979 o066--0.013 

os 0.698 +0.052—0.12 0.8125 o10~0.028 0.8077 5015 0.029 0.812* 0 nos 0.020 
Ho [Km/s/Mpc] 48.8 5 9*1 67.8475 62 1.3 65.853123 67.97 66 1a 
My [eV] « 0.46 « 0.79 « 0.072 « 0.17 « 0.066 « 0.16 « 0.071 « 0.17 

ss l1) re m 

rárag [Mpc] 147.25* $35" 0.64 147.163 5 30-0.55 147.26 * 500.61 147.17 30.0.61 


TABLE III. 68% and 95% CL constraints on various free and derived parameters of the ACDM + XQ + M, scenario for several 
observational datasets. 


Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 
Qeh? 0.1181*5 0930.-0.0060  O-118970,0031.0,0058  O-1187 6.0030" 0.0058 0.1189 5 0032.-0.0058 
Wh? 0.022607 5 05026--0.00049 0-02232+000023—0.00046 0-022490. 9,00024.-0.00047 0-02233" 0.00022 0.00048 

1000mc = 1.04117 00544..0.00086 1-041107 6.00044.0.00086 1.04108 00044 .0.00085 1-041117 00043.-0.00085 
T 0.0484* 5 075.-0.018 0.09467 0.0076 0.015  0-0545t0.0083 0.015  0-055460076—0.015 
ns 0.9703 5 00940.0183  0.9627^0.0088 0.01 0.9069 ' 50091 0.018 0-963175 opas. o.017 
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TABLE IV. 68% and 95% CL constraints on various free and derived parameters of the ACDM + Qx + Neg scenario for 
several observational datasets. 


departures from its reference value, one can probe and 
constrain several extended models both of cosmology and 
particle physics that predict extra dark radiation in the 
early Universe, including the cases of additional neutrino 
species and hot relics beyond the standard model of ele- 
mentary particles [21, 119—129]. 


ways as a result of the tension they have with Planck. 
As concerns the effective number of relativistic degrees 
of freedom, we see that the reference value Neg = 3.044 
is always consistent within one standard deviation for all 
the different data combinations analyzed. Therefore we 
find no significant evidence of deviations from the pre- 
diction of standard model of elementary particles and, in 
any case, in this extended parameter space current cos- 
mological and astrophysical observations still constrain 
additional contributions to ANeg = Neg — 3.044 < 0.4 
(at 95% CL), analogously to the standard case of spa- 
tial flatness.” Actually, as we can see in Figure 3 there 


Table IV summarizes the 68% and 95% CL cosmolog- 
ical constraints, while Figure 3 shows the 68% and 95% 
CL contour plots for different cosmological parameters. 
As for the previous extended model, also in this case the 
CMB data alone suggest a non-flat background geometry 
with the case Qgo = 0 excluded at more than 95% CL. 
Conversely, the inclusion of the other CMB-independent 
datasets eliminates such a preference for a curved space- 
time, giving instead indication for flatness within 1o, al- 


DV] 


5 Notice that while this bound ANeg < 0.4 is tight enough to 
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FIG. 2. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario ACDM + Qx + M, for various observational datasets. 


is no correlation between Neg and Qg. Modifying the 
relativistic energy density in the early Universe leads to 
several implications as it changes the sound horizon at 
recombination which is partly degenerate with the late- 
time geometry. In particular higher values of the effective 
number of relativistic species lead to smaller values of the 
sound horizon at recombination, resulting both in fluc- 
tuations in the amplitude of og and in a preference for 
higher-values of the expansion rate, possibly alleviating 
the Hubble tension. This is not what happens in this ex- 
tended model where the same tensions between early and 


severely constrain several exotic scenarios such as additional neu- 
trino species, it is also too large for probing other extensions of 
the standard model involving extra relativistic species decoupled 
at high temperatures of the order of the top quark annihilation, 
whose contribution ANeg ~ 0.027 is much smaller than the cur- 
rent constraining power. 


late time measurements of Ho discussed in the previous 
cases can be still observed, always with a statistical sig- 
nificance ranging between ~ 50 (CMB only vs SHOES) 
and = 3.30 (when the CMB data are combined with BAO 
and/or Pantheon). Similarly, the different data combina- 
tions provide values of Sg that are always in tension at 
more than 95% CL with each other. 


4. ACDM + Qr + M, + New 


Lastly, we consider an extended non-flat cosmological 
model where we freely vary both the total neutrino mass 
and the effective number of relativistic degrees of free- 
dom. We refer to this model as ACDM + Ox + M, + 
Neg. The results obtained for this model are given in 
Table V at 68% and 95% CL, while the 1D and 2D pos- 
teriors distributions of different cosmological parameters 
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FIG. 3. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario ACDM + Ox + Neg for various observational datasets. 


are shown in Figure 4. Simultaneously varying the to- 
tal neutrino mass and the effective number of relativistic 
degrees of freedom, any significant change is observed 
neither about the Planck preference for a closed Uni- 
verse nor for the Ho-tension. In particular form the CMB 
data alone we get Oxo = —0.074*0029 at 9596 CL and 
Ho = 48.1435 km/s/Mpc at 68% CL, in tension with lo- 
cal measurements at the level of 4.70. As for the previous 
models, including the BAO (Pantheon) likelihood, that 
disagree with Planck in this extended model, the prefer- 
ence for a closed Universe disappears and the constraints 
for the expansion rate become in tension at the level of 
3.30 with the independent result of the SHOES collabo- 
ration. Also in this case the Ho tension drives the values 
of Sg to be in disagreement between 3 and 4 standard de- 
viations with each other, depending on the different data 
combination considered. As concerns the bounds on the 
effective number of relativistic degrees of freedom, also 


in this case we do not find any evidence for a deviation 
from the value expected in the Standard Model and the 
constraints remain basically unchanged with respect to 
the previous case without neutrinos. Similarly, the upper 
bounds on the total neutrino mass are not affected by the 
inclusion of a free effective number of relativistic species 
and they are similar to those obtained within the ACDM 
+ Qg + M, cosmological model. In particular exploiting 
the CMB data alone we can get M, < 0.81 eV at 95% 
CL (with one sigma indication for a total neutrino mass 
different from zero), while including also the BAO and 
Pantheon data this bound can be improved to M, < 0.18 
eV for CMB+BAO(+Pantheon) and M, « 0.16 eV for 
CMB^4-Pantheon, all at 9596 CL. 
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TABLE V. 68% and 95% CL constraints on various free and derived parameters of the ACDM + Qx + M, + Neg scenario 


for several observational datasets. 


B. Non-flat v, CDM. and its extensions 


In this section we describe the extensions of the non- 
flat wCDM model. We mention that for all extensions of 
the nonflat CDM model, we have made use of the same 
datasets and their combinations. Moreover, we also note 
that for each parameter we report either its 68% and 95% 
CL constraints or the upper/lower limits at 68% and 95% 
CL. 


1. wCDM + Ox 


In Table VI we have summarized the constraints on the 
wCDM + Qg model, and the corresponding marginalized 
1D and 2D posterior distributions are shown in Figure 5. 

We notice that for CMB alone, we get an evidence 
of a closed Universe at more than 6896 CL, however, 
within 9596 CL, the spatially flatness of the Universe 
is recovered. The dark energy equation of state is per- 
fectly consistent with a cosmological constant, giving 
w = —1.31*048 at 68% CL for CMB alone. However, 
even if we vary this parameter for CMB alone, unlike in 
phantom dark energy Universe where w < —1 increases 
the expansion rate, here this does not happen and the 
Hubble constant is almost unconstrained, as one can see 
that has a lower mean value Ho = 61159 km/s/Mpc at 
68% CL (CMB alone), but due to the large error bars, the 
tension with the SHOES measurement is reduced down to 
1.2c. This lowering of the mean value is due to the effects 
of curvature that we are witnessing in this case. 

The inclusion of BAO to CMB changes the constraints 
significantly, because of their mutual tension in a curved 
universe, and we see that this scenario becomes very 
close to the spatially flat ACDM one. However, the re- 


sults for CMB+Pantheon are much more interesting. We 
see that for this combination we get a closed Universe 
at more than 95% CL (Qg = —0.028*0020 at 95% CL 
for CMB+Pantheon) together with a phantom Universe 
at more than 95% CL (w = —1.22*011 at 95% CL for 
CMB-+Pantheon), i.e. an indication for a phantom closed 
universe, as already noticed in [11]. Similar to the CMB 
alone case, here too, we get a smaller value of the Hub- 
ble constant (Ho = 61.2734 km/s/Mpc at 68% CL for 
CMB-+Pantheon) and hence a bigger value of Qmo due 
to the positive correlation between Ho and Qmo- 

Finally, if we look at the combination of 
CMB-4BAO-Pantheon (column 5 of Table VI), 
because of the presence of the BAO data that are in 
tension with Planck, we again find that the spatial 
flatness is strongly suggested. At the same time w is in 
agreement within 68% CL with a cosmological constant 
w = —1. There is one point to note as well. Here, 
we see that Ho is mildly increased (Ho = 68.30*034 
km/s/Mpc at 68% CL) than its predicted value within 
the ACDM paradigm by Planck 2018 [8], but the tension 
with SHOES is still at 3.6c. 


2. wCDM + Ox + M, 


Following the same path as the previous section, we 
first consider an expansion which includes the total neu- 
trino mass as an additional free parameter of the sample. 
We refer to this model as wCDM + Ox + M, and we 
summarize the 68% and 95% CL constraints in Table VII, 
showing the correlations among the different parameters 
in Figure 6. 
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FIG. 4. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario ACDM + Ox + M, + Neg for various observational datasets. 


In this case, considering only the CMB data, a spa- 
tially lat background geometry is disfavored at slightly 
more than 9596 CL with the constraint on the curvature 
parameter reading Qo = —0.073*0056 at 95% CL. In- 
terestingly, a strong preference for a curved spacetime is 
found combining CMB and Pantheon (from which we get 
Nko = —0.030*0 025 at 95% CL). Conversely, including 
the BAO measurements, flatness becomes always con- 
sistent within one standard deviation for all the other 
different data combinations, but these results should be 
taken with caution because of the data inconsistency in a 
curved universe. As concerns the dark energy equation of 
state, using only the CMB data we find that both phan- 
tom (w « —1) and quintessential behaviors (w > —1) 
are allowed within the 68% CL results. Combining CMB 
with BAO (+Pantheon) no significant differences are ob- 
served, as well. Conversely, it should be noted that the 
results for CMB+Pantheon (w = —1.281(32 at 95% 


CL) suggest a preference for a phantom dark energy 
at the level of more than 95% CL. As for the baseline 
wCDM + Oy model, also in this extended parameter 
space, allowing phantom models of dark energy does not 
alleviate the tensions, and from CMB+Pantheon thre is 
indication for a phantom closed universe, as in [11]. In- 
deed the value of the Hubble parameter inferred by the 
CMB data, Ho = 527°, km/s/Mpc at 68% CL, and the 
independent measurement provided by the SHOES collab- 
oration are in tension at the level of 3.40. This preference 
for lower values of Ho is mostly driven by the degeneracy 
with the spatial curvature. Indeed, when the BAO data 
are combined with the CMB, no evidence is found for cur- 
vature and the constraints on Ho look very different and 
become similar to the usual results obtained within the 
flat ACDM model. Conversely, for CMB+Pantheon, due 
to the preference for a curved spacetime of this dataset, 
we still obtain a lower value Ho = 60.2 + 2.6 km/s/Mpc 
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TABLE VI. 68% and 9596 CL constraints on various free and derived parameters of the wCDM + Qx scenario using several 


observational datasets. 


at 68% CL, in strong tension with local measurements at 
about 4.6c. Finally, it is worth pointing out that in this 
model the bounds on the total neutrino mass obtained 
using the Planck measurements (M, « 0.83 eV at 9596 
CL) remain almost unchanged with respected to the pre- 
vious scenarios, while the bounds for CMB+Pantheon 
are significantly relaxed to M, « 0.46 eV at 9596 CL. 
Anyway, the most significant impact on the bounds on 
relic neutrinos is obtained by combining CMB and BAO 
from which we obtain M, « 0.19 eV at 9596 CL. Con- 
sidering CMB--BAO-rPantheon all together this bound 
is robust and almost unchanged. 


3. wCDM + Ox + Neg 


Now we analyze the CDM + Os + Neg model, vary- 
ing the effective number of relativistic degrees of free- 
dom as additional parameter instead of the total neutrino 
mass. The results for this model are given in Table VIII, 
and the triangular plots are shown in Figure 7. 

Notice that replacing the total neutrino mass with the 
effective number of neutrino species the results for the 
curvature parameter are changed and the Planck pref- 
erence for a closed Universe is reduced. In this case 
the Planck limit on curvature Oxo = —0.049*0093 at 
95% CL is consistent with flatness at the level of 1.1c. 
Interestingly, CMB+Pantheon combined together pre- 
fer instead a closed Universe and we obtain (Qo = 
—0.028*0020 at 9596 CL. Such a preference disappears 
once we introduce the BAO measurements, as noticed 
before. As concerns the Ho-tension, we see that, due to 
the correlation between the expansion rate and Neg the 
CMB bounds are less tight, but they show in any case a 


preference for smaller values Ho = 593, km/s/Mpc at 
6896 CL. Due to the large uncertainty, the tension with 
the SHOES estimation below 2c. Combining CMB+BAO 
instead we get Ho = 68.2* 15 km/s/Mpc at 68% CL, close 
to the value obtained within the flat ACDM model but 
reducing the tension with SHOES at the level of ~ 2.20, 
due to the larger error-bars. On the other hand, consid- 
ering CMB+Pantheon, the preference for a curved Uni- 
verse lowers the value of Hg to Hp = BLl km/s/Mpc 
(at 68% CL) while the uncertainty remains almost un- 
changed. As a result, the tension reaches the level of 
4.80. Notice also that these tensions strongly affect the 
bounds on the matter-density parameter Qmo and con- 
sequently the value of Sg since both of them are pushed 
toward higher (lower) values when lower (higher) Ho are 
preferred by the different data combinations. Concern- 
ing the equation of state of Dark Energy, the CMB mea- 
surements give a very relaxed bound on w that is com- 
pletely consistent with a cosmological constant. This 
result remains true also combining the CMB and BAO 
data, while, as in the previous model, the combination 
CMB-+Pantheon strongly suggest phantom dark energy, 
excluding a cosmological constant at more than 95% CL. 
Finally, it should be noted that the bounds on the ef- 
fective number of relativistic degrees of freedom remain 
basically unchanged and, exactly as in the case of (flat 
and curved) models with a cosmological constant, addi- 
tional contribution are constrained to ANeg < 0.4 for all 
the different datasets. 
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FIG. 5. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario wCDM + Qx for various observational datasets. 


4. wCDM + Ox + M, + Neg 


In this last extension of wCDM + Qg, we simulta- 
neously vary the total neutrino mass and the effective 
number of relativistic degrees of freedom. We refer to 
this model as wCDM + Ox + M, + Neg providing the 
results in Table IX, and showing the 1D and 2D poste- 
rior distributions of different cosmological parameters in 
Figure 8. 


In this extended scenario the Planck preference for 
a closed Universe is reduced below 2e and the CMB 
data constrain Ngo = —0.067*0005 at 9596 CL. Con- 
versely, a preference for a curved geometry is found an- 
alyzing CMB+Pantheon, as in this case we get Qgo = 
—0.029*0029 at 9596 CL. This is the same behavior ob- 
served studying the previous extensions. Indeed also the 
results for Ho remain stable with respect to the case 


without neutrinos. In particular from the CMB mea- 
surements we get Ho = 54*2» km/s/Mpc at 6896 CL, 
reducing the tension with the SHOES collaboration value 
at less than ^ 20 due to the large error, and not to 
a better overlapping. Combining the CMB data with 
BAO or Pantheon we obtain Ho = 68.2*55 km/s/Mpc 
and Ho = 60.3*25 km/s/Mpc both at 68% CL, respec- 
tively, producing basically the same tensions among the 
different datasets discussed in the previous model. As 
concerns the equation of state of the dark energy, ana- 
lyzing CMB+Pantheon we still observe a preference at 
more then 9596 CL for phantom models, while all the 
other datasets are in agreement with the cosmological 
constant within one standard deviation. As concerns the 
constraints on the total neutrino mass, the upper bounds 
are not changed significantly and the most constraining 
bound obtained for CMB+BAO (M, « 0.20 eV at 95% 
CL) is only slightly relaxed varying also Neg in the sam- 
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TABLE VII. 68% and 9596 CL constraints on various free and derived parameters of the wCDM + Ox 4 


several observational datasets. 
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TABLE VIII. 68% and 95% CL constraints on various free and derived parameters of the wCDM + Ox + Neg scenario using 
several observational datasets. 


ple. Similarly, no significant change is observed in the 


observational constraints for Neg. 


C. Non-flat wow, CDM. and its extensions 


Here we discuss the observational constraints on the 
non-flat t 9w4 CDM model and its various extensions 
through the inclusion of the neutrino sector. As writ- 
ten before, we have used exactly the same datasets and 


their combinations and for each parameter we report ei- 


ther its 6896 and 9596 CL constraints or the upper/lower 
limits at 68% and 95% CL. 


1. wows4CDM + lk 


In Table X we have summarized the observational con- 


straints on the free and derived parameters of this model 


and in Figure 9 we show the 1D and 2D posterior distri- 
butions on a few interesting parameters. Let us describe 
the impacts of CMB alone and its combination with other 
non-CMB datasets. 
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FIG. 6. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario wCDM + Qx + M, for various observational datasets. 


As displayed in the second column of Table X, we see 
that CMB alone indicates the preference of a closed Uni- 
verse at more than 68% CL, but Qx is in agreement with 
zero within 2e (Qg = —0.038*0 p75 at 68% CL). And as 
usual we have a lower mean value of Ho for this case 
(Ho = 63*19 km/s/Mpc at 68% CL) but with very large 
error bars. Due to such high error bars, the tension with 
SHOES is reduced down to ~ lo. The current value of 
the dark energy equation of state even though has its 
mean value in the phantom state (wọ = —1.3*,2 at 68% 
CL) within 68% CL is in agreement with the cosmological 
constant. Also, the dynamical character of the dark en- 
ergy equation of state quantified through wa (= —0.6*12 
at 68% CL) is in agreement with zero, i.e. the dynamical 
DE is not favoured. 

When BAO data are added to the CMB, the possi- 


bility of nonzero curvature of the Universe is present 
only at 68% CL (Qk = —0.0042*00026 at 6896 CL for 


CMB+BAO) and the usual shift of Ho towards lower 
mean values is present as well, with a reduction of its 
error bars with respect to the CMB only case (Ho — 
64.615 km/s/Mpc at 68% CL), restoring the tension 
with local measurement at ~ 3.90. Concerning the free 
parameters wo and Wa we see that wo is strictly in the 
quintessence regime at more than 9596 CL while wa is 
non-zero at more than 9596 CL, that means we have an 
evidence of the dynamical quintessence dark energy for 
this case at more than 2c. 


When Pantheon data are added to CMB, we recover 
the evidence of nonzero curvature at more than 9596 CL 
(Qe = —0.031*002? at 9596 CL), while the tension of 
Ho with SHOES exacerbates. About the properties of 
the dark energy parameters, we see that the mean value 
of wo is perfectly in agreement with the cosmological 
constant wo = —1 ( wo = —1.04*016 at 68% CL for 
CMB-+Pantheon), but the dynamical nature of the DE is 
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FIG. 7. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario wCDM + Qx + Neg for various observational datasets. 


signaled through Wa which is found to be nonzero within 
68% CL (wa = —1.2*9? at 68% CL for CMB+Pantheon) 
while it becomes consistent to zero at 95% CL. 

For the combined analysis CMB+BAO+Pantheon we 
are able to break all of the degeneracies and have a very 
mild improvement of the Hubble constant agreement, 
ie. Ho = 68.0149'32 km/s/Mpc at 68% CL for this 
case. The curvature parameter is consistent with zero at 
slightly more than 68% CL. Additionally, we recover the 
quintessence nature of wo at slightly more than 68% CL, 
similar to what we observed for the CMB+BAO analy- 
sis, and the same we can say for wa, that is only found 
to be nonzero at more than 6896 CL, showing a hint for 
a dynamical nature, while within the 9596 CL wa = 0 is 
back in agreement with the data. 


2. wowa CDM + OK + My, 


As usual, we start extending the parameter space of 
the model considering massive neutrinos. We refer to 
this model as wow, CDM + Ox; + M, and we show the 
cosmological constraints on the parameters in Table XI 
and in Figure 10. 

Considering only the CMB data the value found for the 
curvature parameter is Qgo = —0.065*008$ at 95% CL 
and one can see that while flatness is disfavored by the 
6896 CL limit it is consistent within the 9596 CL bounds, 
similarly to the baseline case discussed above. Interest- 
ingly, including BAO measurements, we still find a pref- 
erence for a closed Universe at the level of 6896 CL even 
though flatness is still consistent with the 9596 CL lim- 
its. Conversely, when the Pantheon data are added to the 
CMB measurements, we recover the evidence of nonzero 
curvature at more than 95% CL (Qo = —0.034*0 054 at 
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Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 
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mo TTE. E e n 

os 0.76* 517 0.22 0.818 5 25 "o. 054 0.811* 0/020. 0.058 0.814* 016 0.037 

Ho [Km/s/Mpc] 54* 12723 68.252 58 60.3* 5 9*5 67.8*12* 3 

My [eV] « 0.49 « 0.81 « 0.084 « 0.20 « 0.203 « 0.41 « 0.080 « 0.20 
Nest 3.0475 200.39 2.94 5 900.37 3.02 5 25 0.37 2.9479 190.39 
Ss 0.984* 0071 -0.16 0.82375 016.-0.035 0.927 ^ 0038 - 0.076 0.82315 015.-0.023 
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TABLE IX. 6896 and 95% CL constraints on various free and derived parameters of the wCDM + Nx + My + Neg scenario 
using several observational datasets. 


Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+ Pantheon 
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Nko —0.038* $015. 0.052 —0.0042* $536. 0.0066 —0.031*$ oii ^ 0.021  —0.0032* 5 oo30.-0.0065 
Qmo 0.417524" 0.30 0.34275 022" 0.047 0.382* 5035 0.058 0.3081000760015 
» UT M c MEETS UL NETT S 
Ho [Km/s/Mpc] 63* 15 24 64.6 y 5 pa 61.0 $$ 51 68.011035 1.7 
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TABLE X. 68% and 9596 CL constraints on various free and derived parameters of the wow; CDM + Qx scenario using several 


observational datasets. 


95% CL for CMB+Pantheon). As concerns the Hubble 
parameter, from the CMB data we obtain a lower mean 
value with respect to the baseline case without neutri- 
nos and our 95% CL now reads Ho = 55*21 km/s/Mpc. 
Due to the large error-bars this result is in agreement 
with SHOES within 2c. On the other hand, combining 
CMB+BAO we find Hy = 64.6712 km/s/Mpc at 68% 
CL and in this case, due to the smaller uncertainty, the 
tension is increased to 3.90. Finally, when we combine 
CMB-+Pantheon, the combined effects of curvature and 
neutrinos, lead to smaller values for the expansion rate, 
Ho = 59.7734 km/s/Mpc, at 68% CL and the result is 


in strong tension with local estimation at 5c. The dif- 
ferent datasets analyzed in this extended model provide 
discordant predictions also for the dark energy sector. In 
particular, exploiting only the Planck satellite measure- 
ments wo turns out to be unconstrained while only an 
uninformative upper bound wa < 1.9 at 9596 CL can be 
derived for the parameter who carries information about 
the dynamical nature of DE. Conversely, including also 
the BAO measurements we find both a strong prefer- 
ence for a quintessential dark energy (wo — —0.51*048 
at 9596) with a dynamical behavior, as the upper bound 
Wa < —0.4 rules out the case wa = 0 at 95% CL. On 
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FIG. 8. 1-dimensional marginalized posterior distributions and the 2-dimensional joint contours for the most relevant parameters 
of cosmological scenario wCDM + Ox + My + Neg for various observational datasets. 


the other hand, CMB+Pantheon is again in agreement 
with a cosmological constant within 68% CL, and the 
upper bound wa < 0.24 at 95% is in agreement with a 
non-dynamical equation of state for the dark energy. Fi- 
nally, combining CMB+BAO+Phantom all together we 
find that quintessential dynamical models are preferred 
at more than lo but both phantom models and non- 
dynamical dark energy are allowed within the 95% CL 
constraints. As concerns neutrinos, the limits on their 
total mass are not drastically changed by the modifica- 
tion in the dark energy sector and the 95% CL bound ob- 
tain for CMB+BAO (M, < 0.22 eV at 95% CL) is only 
slightly larger than the respective limit obtained within 
curved ACDM or wCDM cosmologies. This full dataset 
is completely in agreement with a flat universe. 


3. WoWa CDM + Ox + Neg 


We now allow the effective number of relativistic de- 
grees of freedom to freely vary in curved cosmologies with 
dynamical models of dark energy. We refer to this sce- 
nario as wow; CDM + Ox + Neg. The 68% and 95% CL 
constraints on the parameters of the model are given in 
Table XII and in Figure 11 we show the marginalized 1D 
and 2D posterior distributions of the parameters. 


Considering the CMB bounds, also in this case we find 
that flatness, while disfavored by the 68% CL limit, is 
consistent within the 95% CL bounds. Instead, for the 
Hubble constant we get Ho = 6415) km/s/Mpc at 68% 
CL that, due to the large errors, is consistent with the 
SHOES measurement within one standard deviation. No- 
tice also that in this model the parameters involving mod- 
ifications to the dark energy sector are basically uncon- 
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FIG. 9. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by wowaCDM + €; using various observational datasets. 


strained by the CMB measurements. Including the BAO 
data a curved background geometry is always slightly 
preferred at 6896 CL, but flatness is consistent within 
the 9596 CL limits, as well. The value of the Hubble con- 
stant Ho = 64.6* 21 km/s/Mpc at 68% CL remains stable 
with respect to the CMB result, but with a much smaller 
uncertainty. Consequently, the tension with SHOES is 
increased at the level of ~ 3.60. Concerning the dark en- 
ergy parameters, for this dataset quintessential models 
are preferred and the constraints on wọ reads at 9596 CL 
wo = —0.55*0 57. Instead the upper bound w, « —0.1 
excludes a non-dynamical evolution at 9596 CL. Com- 
bining CMB and Pantheon we find the usual evidence 
for a curved Universe (Qj = —0.031 + 0.021 at 95% 
CL) resulting into a preference for smaller values of the 
expansion rate (Ho = 60.9*2:3 km/s/Mpc at 68% CL). 
This is in large tension with SHOES at more than 4.60. 
Furthermore in this case both phantom and quintessen- 


tial models of dark energy are allowed, or in other words 
wo is in agreement with a cosmological constant, and 
the bounds on Wwa do not give substantial information 
about its dynamical evolution. When instead we com- 
bine CMB+BAO-+Pantheon all together, we do not find 
evidence for curvature and for then expansion rate we 
obtain Ho = 67.6 + 1.3 km/s/Mpc at 68% CL. This 
result is similar to the best fit value that can be ob- 
tained within the standard flat ACDM model. For the 
same dataset quintessential dynamical models of dark en- 
ergy turns out to be preferred at 68% CL even though 
both phantom models and non-dynamical evolution are 
allowed at 95% CL and the bounds on the dark energy 
parameters reads wg = —0.89*0:22 and wa = —0.75 + 1.1 
at 9596 CL, respectively. Finally, concerning the con- 
straints on the effective number or relativistic degrees of 
freedom, it is worth pointing out that they remain al- 
most unchanged with respect to the other curved models 
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Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 
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Ho [Km/s/Mpc] 55 i21 64.6 56 23 59.7 $ 008.2 67.98 851.7 


M, [eV] 0.42653 < 0.84 < 0.092 < 0.22 < 0.272 < 0.49 < 0.090 < 0.21 
+0.090+0.14 +0.022+0.038 +0.036+0.070 +0.018+0.036 
Sg 0.977 9.067 0.16 0.854" 9017 0.041 0.938" 5 956. 0.075 0.8407 9018. 0.036 
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TABLE XI. 68% and 95% CL constraints on various free and derived parameters of the wowa CDM + Ox + Mz scenario using 


several observational datasets. 


of dark energy analyzed in this work for all the different 
data combinations. 


4. WoWa CDM + Or + M, + Neg 


In this last extension we simultaneously vary the total 
neutrino mass and the effective number of relativistic de- 
grees of freedom. We refer to this scenario as wowa CDM 
+ Ok + M, + Neg, providing the results in Table XIII, 
and the contour plots in Figure 12. 

In the large parameter-space allowed in this model, 
the results for curvature does not change significantly. 
In particular, as in the previous extensions, the CMB 
data always prefer curved spaces at 6896 CL, while flat- 
ness is always allowed within the 95% CL results (Q xo = 
0.056*0025). As concerns the Hubble parameter, CMB 
data give Hg = 57*19 km/s/Mpc at 68% which is sim- 
ilar to the bound obtained without including the effec- 
tive number of relativistic degrees of freedom in the sam- 
ple. Indeed, as we already pointed out, the preference for 
smaller values of Ho is mostly driven by the combined ef- 
fects of curvature and massive neutrinos and adding Neg 
does not reduce the tension between the different datasets 
discussed so far. Concerning dark energy, both wo and wa 
turns out to be basically unconstrained for these dataset 
and we can derive only a 9596 CL upper bound wa < 1.6 
which is quite uninformative about the possible dynam- 
ical nature of the late time expansion. Including BAO 
in the picture, we basically recover the same results dis- 
cussed without varying Neg, as well. In particular we 


find Qo = —0.0037*00072 at 9596 CL and it should be 
noted that, while this bound is very shrink around zero, 
flatness is still slightly disfavored at 68% CL. Similarly 
for Ho we find 64.4*21 at 68% CL producing the same 
tensions discussed in the case without Neg. On the other 
hand, for the constraints wo = —0.52* 01s and wa < —0.4 
(both at 95% CL) we can observe a strong indication for 
quintessential models of dark energy with a dynamical 
evolution that are indeed preferred at more than 9596 
CL. When CMB+Pantheon are combined together, we 
recover the usual evidence for a closed Universe at more 
than 95% CL (Qo = —0.033+21 at 95% CL) and we get 
a lower value Ho = 59.8*25 km/s/Mpc at 68% CL for 
the Hubble parameter, in strong tension with local mea- 
surements as usual. In this case both quintessential and 
phantom models of dark energy are both allowed within 
the 68% CL (wo = —1.09* 9:2. at 68% CL), i.e. wo is con- 
sistent with a cosmological constant, and only an upper 
bound wa < 0.24 can be derived on wa at 95% CL. On 
the other hand, combining CMB--BAO Pantheon all to- 
gether the evidence for a closed Universe disappears, and 
the value of the expansion rate is shifted towards higher 
values HO = 67.6 + 1.2 km/s/Mpc at 68% CL, similar 
to the standard ACDM case. Furthermore in this case 
wo is consistent with a cosmological constant within the 
68% CL, and wa = 0 within the 95% CL. Finally, con- 
cerning the mass or relic neutrinos, exploiting the CMB 
temperature and polarization measurements only, we get 
M, « 0.753 eV at 9596 CL, while combining the CMB 
and BAO data we can improve this limit to M, « 0.231 
eV at 9596 CL. This bound remains almost unchanged 
including also Pantheon data and the results are very 
similar to those discussed for the previous cosmological 
models. The same discussion applies also to the effec- 
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FIG. 10. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by wow; CDM + Qx + My using various observational datasets. 


tive number of relativistic degrees of freedom and in par- 
ticular we can appreciate that additional contribution 
always constrain AN.g < 0.4, independently from the 
parametrization introduced in the dark energy sector. 


D. Non-flat PEDE and its extensions 


Here we describe the final model of this article and 
its extensions. In a similar fashion we have considered 
exactly the same datasets and their combinations and 
concerning the observational constraints, we report either 
their 6896 and 9596 CL constraints or the upper/lower 
limits at 6896 and 9596 CL. 


1. PEDE + Qk 


In Table XIV we have summarized the observational 
constraints on this scenario, and in Figure 13 we show the 
1-dimensional posterior distributions of some parameters 
and the 2-dimensional joint contours between several pa- 
rameters. 


We have some interesting results in this case. We see 
that all the datasets indicate the preference of a nonzero 
spatial curvature. Precisely, we observe that Qg remains 
nonzero at more than 95% CL for CMB, CMB+ Pantheon 
and CMB--BAO--Pantheon while for CMB+BAO, Ox 
remains nonzero at more than 6896 CL. Additionally, for 
all the cases explored here, the mean values of Qg are 
negative which indicate the preference of a closed Uni- 
verse. 


Concerning the estimated values of Ho, we see that 
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TABLE XII. 68% and 95% CL constraints on various free and derived parameters of the wow; CDM + Ng + Neg scenario 
using several observational datasets. 


Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 
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Qo 0.05675 oio o.oro  —0-0037 0 0035.-0.0008  —0-093-0.011 0.021  —0-0025 f Qo33 0.0005 
Qmo 0.5170 3170.39 0.342 9980-047 0.403003940071 930999052" 0.016 
^ IINE UT MEUS -o 
Ho [Km/s/Mpc] 57 1573s 644757 4.6 59.8 55 55 07.012121 
M, [eV] « 0.365 « 0.753 « 0.092 « 0.231 « 0.279 « 0.49 « 0.098 « 0.22 
Nar UE. 29001000 path oe gost) uo 
s. Tu MEE UI MEET CU EMEND. 
raras [Mpc] 147.397 55 5 76 147.94 555 26 147.577 938 148.277 93.6 


TABLE XIII. 68% and 95% CL constraints on various free and derived parameters of the wow; CDM + Qx + My + Neg 
scenario using several observational datasets. 


CMB alone predicts a very low value (Ho 
km/s/Mpc at 68% CL) and hence the tension is about 
~ 3.50. However, the inclusion of BAO quite surpris- 
ingly increases the mean value of Ho and decreases the 
error bars leading to Hp = 70.57 ^94; km/s/Mpc at 68% 
CL and thus the tension with SHOES is reduced down 
In the contrary, the inclusion of Pantheon to 


CMB slightly increases the Hubble constant compared 


to 2c. 


+0.73 


5813$ 


to CMB alone leading to Hp = 62.3 + 2 km/s/Mpc at 
68% CL (CMB+Pantheon) and again the tension with 
SHOES is found to be at ~ 4.80. So, effectively for CMB 
and CMB+Pantheon we can see a clear tension. How- 
ever, similar to the previous case with CMB--BAO, the 
combined analysis CMB--BAO-4-Pantheon also leads to a 
slightly higher value of Ho ( = 69.94* 079 
68% CL for CMB+BAO-+ Pantheon) compared to Planck 


km/s/Mpc at 
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FIG. 11. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by wowaCDM + Qg + Neg using various observational datasets. 


alone (within ACDM) and the tension with SHOES is re- 
duced down to ~ 2.50. 


2. PEDE + Ox + M, 


As usual, we start expanding the parameter space in- 
troducing massive neutrinos in the cosmological model. 
We summarize the results obtained in this case in Ta- 
ble XV and in Figure 14. 

Including neutrinos we find that the CMB and 
CMB+Pantheon data always prefer a closed Universe 
at more than 95% CL, while this preference is re- 
duced to 68% CL both for CMB+BAO and for 
CMB+BAO-+Pantheon. As concerns the Hubble param- 
eter, in this case, due to the degeneracy with massive 


neutrinos, the CMB data alone show a preference for 
lower values giving Ho = 51.2'97 km/s/Mpc at 6896 
CL. Despite the large error-bars, this value is in tension 
with SHOES at 3.20. Interestingly, including BAO mea- 
surements, we find back the same preference for larger 
values discussed in the baseline case with the 6896 CL 
result reading Ho = 7010-45 km/s/Mpc. In this case the 
tension with the SHOES result is reduced down to 2.50. 
Conversely for CMB4-Pantheon, we get Hp = 62.3 + 1.9 
km/s/Mpc at 68% CL and the tension with SHOES is 
back to be statistically significant. Finally combining 
CMB+BAO-+Pantheon all together we obtain an inter- 
mediate value Ho = 69.78191} km/s/Mpc at 68% CL 
which is very similar to the result that one can derive 
within the standard ACDM model of cosmology. Notice 
that, as already discussed for all the previous models, 
also in this case the Ho tension has an impact on the 
estimation of the matter content of the universe, result- 
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FIG. 12. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by wowa CDM + Ox + M, + Neg using various observational datasets. 


ing into a tension for Sg. As concerns relics neutrinos, 
the constraints on their total mass remain almost un- 
changed with respect to the previous scenarios analyzed 
in this work, but it is worth noting that in this case 
we observe an improvement in the constraining power 
of CMB+Pantheon that give M, < 0.26eV at 95% CL. 
When CMB and BAO are combined together we can im- 
prove this upper limit to M, < 0.22 eV at 95% CL, while 
combining CMB-4-BAO-4 Pantheon all together we get a 
less tight bound M, « 0.37 eV at 9596 CL. 


3. PEDE + Ox + Neg 


We now replace massive neutrinos with the effective 
number of relativistic degrees of freedom Neg in the cos- 


mological model. We summarize the results obtained 
within this scenario in Table XVI and in Figure 15. 


Analyzing the Planck satellite measurements of CMB 
temperature anisotropies and polarization only, in this 
case we observe an indication for a closed Universe at 
more than 9596 CL, with the curvature parameter read- 
ing Oxo = —0.038*0025 (at 9596 CL). Such an indica- 
tion remains true for CMB+Pantheon while combining 
CMB--BAO and CMB+BAO+Pantheon it is reduced 
to the level of 68% CL, with flatness always allowed 
within the 9596 CL results. Similarly to the baseline case, 
from the CMB we get a 68% CL value Ho = 57.9 10 
km/s/Mpc for the Hubble parameter which is about 3.70 
in tension with the independent local measurement pro- 
vided by the SHOES collaboration. Interestingly, also 
in this case, combining CMB+BAO the tension is re- 
duced at 1.90 due to the higher value Hp = 70.0 + 1.2 
km/s/Mpc at 68% CL preferred for this dataset. Con- 
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Parameters CMB CMB+BAO CMB-+ Pantheon CMB+BAO-+ Pantheon 
QA OLLBITOOBETU N QoS torao GUN —qaigat9.0OlFO 0088 — 0,1203F0:001770-0078 
Mh? —— 0,0226078,00017+0-00082 0 qpn4n0 00018 +9-00020 0 oang5+0-00010+0.000828 0 gpny7+0.00018 +0-000%0 
1000mc 1.041164 5 0,035 0.00066 1-040987 540031 0.00060 1.04112 5, 00032-.0.00064 1040904000031- 0.00062 

r UID MET TUI PEE 
ns 0.97067 0047--0.0094 0.9662 TG 0044-0.0084 0.9095 * 5 0045..0.0089 0-906487 o 0044.-0.0090 

VEU I EE D EEUU EU 
- Eur MEET UID META ur  —0.0038+0:0020 0.0034 
Qmo 0.42675 0690.12 0.2866" 5 0065. 0.013 0.366 5 024 0.042 0.2930* 50062. 0.012 
as 0.8107 $019. 0.038 0.852 5 o084 0.017 0.826* $611 0.021 0.853075 0840.07 

Ho (Km/s/Mpe] 58.0435154 T ean i 69.94 +0: OH 
Ss 0.961* 0 o9. 0.096 0.833 ty o3. 0.025 0.912* 0 025 0.043 0.843^ 5/912 0.024 
rárag [Mpc] 147.34 * y 35 0.59 147.14* 5/51 0.60 147.308 30-058 147.03* 5:35 5:60 


TABLE XIV. 68% and 95% CL constraints on various free and derived parameters of the PEDE + Qx scenario using several 
observational datasets. 


Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 

Qh? a R e EXIT EC EET RES: 
"EET UD Xu E qu uu 
I000.c 1041007 a ioigggf? 990889 90002  pirpat 9908249.00002 ] pner 90088-+0.00061 

: Xu MET un MET n MEE 

ECT n MET ETT UEmM T. 

M MEME EET MEE ET 
ma -0084*9 04i o e  _o.on15+0.0068+0.011 o qnag+0-0028+0.0008 

so II MEME T MEM du EP eI. 

! ^ LE Gu MEME MEE UU MEET S 
Moo Ho [km/s/Mpq 51.2873 DET. 92.319527 oann 
M, [eV] < 0.459 < 0.81 < 0.095 < 0.22 < 0.123 < 0.26 < 0.167 < 0.37 

s. gosatoosejon.  ogagtomrioosi  poogtgoastenso paaytoiae roosa 

raag [Mpe] 147.2503203 qay agtositdeo ——— 1472915287938  qqy.g1+03040.62 


TABLE XV. 68% and 95% CL constraints on various free and derived parameters of the PEDE + Ox + M, scenario using 
several observational datasets. 


versely CMB+Pantheon prefer again smaller values of same that seems can solve the “CMB tension” between 
the Hubble parameter (Ho = 62.3*1* km/s/Mpc at 68% the different CMB datasets [130]. 

CL) in contrast with SHOES at 5.lo. Finally combin- 

ing CMB+BAO+Pantheon all together we find an inter- 

mediate value between the two cases which, being very 

similar to the value obtained within the standard cos- 

mological model, is in tension with local measurements 


as well. As concerns the effective number of relativistic 4. PEDE + Qr + M, + Nen 

degrees of freedom, we see that also in this case the con- 

straints are almost unchanged with respect to the results Finally, we simultaneously vary the total neutrino mass 
derived analyzing different parametrizations of the dark and the effective number of relativistic degrees of free- 
energy sector and additional contributions to dark radi- dom, reporting the results in Table XVII and in Fig- 
ation are always constrained to be ANeg S; 0.4 at 95% ure 16. 

CL. Anyway it is worth pointing out that for the data Exploiting only the Planck observations of the CMB 


combination CMB+BAO-+Pantheon, the reference value radiation, in the largest parameter space sampled in this 
Neg = 3.044 is disfavored at 68% CL and a preference model, we obtain a preference for a curved cosmological 
for smaller values is observed. This smaller value is the spacetime at more than 9596 CL, obtaining a very small 
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FIG. 13. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by PEDE + Qx using various observational datasets. 


value of the expansion rate Ho = 50.3*65 km/s/Mpc 
at 9596 CL. The preference for a such small value is al- 
ways due to the combined effect of curvature and the 
degeneracy with the mass of relic neutrinos that pro- 
duce also a large uncertainty. Combining the CMB and 
BAO data, the preference for a closed Universe is re- 
duced at slightly more than 6896 CL and flatness be- 
comes consistent at 9596 CL. On the other hand the 
constraints on Ho are drastically changed and, like for 
the previous extensions, we instead obtain preference for 
higher values (Ho = 69.9 + 1.3 km/s/Mpc at 68% CL) 
reducing the tension with the SHOES measurement at 
^ 1.90. Conversely, combining CMB and Pantheon we 
find back a strong preference for a curved spacetime (with 
the curvature parameter reading Ngo = —0.021*0 015 at 
9596 CL) and smaller values of the Hubble parameter 
(Ho = 62.4*1* km/s/Mpc at 68% CL) increasing again 
the tension between the different estimations of Hy and 


Ss basically in the same way, and within the same sta- 
tistical significance reported for the previous models. In- 
terestingly, in contrast with the other extensions of this 
baseline model, here combining CMB-4-BAO-4-Pantheon 
all together the preference for a closed Universe disap- 
pears and flatness is always consistent within the 68% 
bound. On the other hand for the Hubble parameter we 
derive results similar to those obtained within the stan- 
dard ACDM model, namely Ho = 68.9715 km/s/Mpc 
at 6896 CL. Notice that also in this extended scenarios 
large disagreements among the different datasets are ob- 
served and adding additional parameters to the baseline 
case does not lead to any significant reduction of such 
tensions. Finally, concerning the neutrino sector, also in 
this case the CMB bound M, « 0.84 eV at 9596 can be 
improved up to M, « 0.26 eV and M, « 0.24 eV in- 
cluding the Pantheon and BAO measurements, respec- 
tively. Conversely, combining CMB+BAO-+ Pantheon 
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FIG. 14. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by PEDE + Ox + M, using various observational datasets. 


the resulting bound is less constrained (M, < 0.35 eV). 
These results are almost the same derived without vary- 
ing the effective number of relativistic degrees of freedom 
and, similarly, also the bound on Neg are not affected 
by the inclusion of neutrinos. In particular, combining 
CMB+BAO+Pantheon we find the same slight prefer- 
ence for smaller values of Neg with respect to the pre- 
diction of the standard model of elementary particles. 
Anyway everything remains consistent with the standard 
model within the 95% CL as well as additional contribu- 
tions are always constrained to be ANeg < 0.4 at 95% 
CL for all the different dataset. 


V. SUMMARY AND CONCLUSIONS 


Modern cosmology has witnessed a remarkable success 
over the last several years and the credit goes a large 
amount of observational data coming from different as- 
tronomical surveys. With the increasing sensitivity in 
the astronomical data, undoubtedly, we have been able 
to put stringent constraints on the cosmological parame- 
ters, and at the same time we are facing new challenges. 
The Ho and Sg tensions have emerged as two biggest 
challenges for the standard ACDM cosmology demand- 
ing its revision, and the Planck preference of a closed 
universe model (see Refs. [8-11, 15, 16]) further argues 
that the curvature of the universe can play a vital role 
in the estimations of the cosmological parameters. Even 
though the earlier investigations were in agreement with 
a spatially flat universe, however, there is no reason to 
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Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+Pantheon 
Meha —— OAITATEOOBBTO OOS  0.1182+0:0030+0 000E 0.118240003040 0083 — 0,1176F9:0032 F0 0000 
E a E a A E aoneIrS opi 

1008aro  1.04120+0:00042+0:00088 1 04115+0:00044+0:00090 1041 16+0:00043+0:00088 1,1429 0.00044 0.00006 
ro n oe a 005105 09001 
ns 0.970275 o090~0.0177 0-9624} 0.0090 0.0181  0-9686t0.0088—0.018 0.95737 0.0097.-0.017 

maota,  aozstomarooe MEET gogatootrenoas 3,oazTo MT. 
Qo —0.038* 5 130.032. —0.0029^ 0 5020.-0.0038 —0-021810:00530.012  — 0.0029 0, Q020.- 0.0039 
Omo 0.428 ^ 069 -0.12 0.2883 0 Qo69 0.014 0.3657 0:021 0.040 0.2956 0 oo69.-0.013 
rs Lu UMEN M s 
Ho [Km/s/Mpc] 57.9 768.0 70.0*12*28 62.3 5 1 $5 68.9* 1 3*5 
Nest 3.0375 19-0 38 2.9479 5070.40 3.0275 940.38 2.85 5 2170.38 
s: lo HE E o 
Tarag [Mpc] 147.567 185 3.66 148.17* 94 4.02 147.5* 197 $7 149.0* 1$ $5 


TABLE XVI. 68% and 95% CL constraints on various free and derived parameters of the PEDE + Ox + Neg scenario using 
several observational datasets. 


Parameters CMB CMB+BAO CMB-+Pantheon CMB+BAO-+ Pantheon 
Neh? 0.1182t0.00s0—0.00s9  0-1181Ż0:0031+0.0060  0-118210.0030—0.0059 0-117810.0029—0.0058 
wh? 0.022547 0.00026—0.00053 0-02232Ż0.00023—0.00044 0-022517 8.00025-0.000450 0:0221810.00024— 0.00045 

1000mc  1.04102'5 00045..0.00092 1.0411510:00045.-0.00086 1:0411210.00047—0.00082 1:0411610.00046—0.00091 
T 0.0476Ż0.0080—0.018 9053575 Sors-o.016 0-051075 0075 0.016 0-0525 10.0081 -0.015 
ns 0.9693t0.0096—0.0196 0.9618 5 Sose—o.01rr 0-9079 0,0094 .0.018 0-9568 10.0086—0.017 

In(10 As) 3.025 18 - 0.042 3.038 0 015. -0.036 3.033" o1 0.037 3.035 0 020 0.038 
Oxo 0.06975 555 -0.ora  —0.0025750025 0.0046  —0.020678.0061.0.013  —0-0017 "9 0032 "0.0084 
Du 0.6170111038 gq. agqo#2-0074+0.018 ggg 5 #0.021+0-043 ag. pggg¥0.007840.017 

os apen a og WI ou enn METH 
Ho Km/s/Mpc| — 50355551 69.9*1 3425 62.475159 68.9 11H24 
M, [eV] 0.89*0-12 < 0.84 « 0.103 « 0.24 « 0.128 « 0.28 « 0.172 « 0.35 
Nett 3.051520. 0.40 2.95 19 0.38 3.021020. 0.30 2.8679 300.38 
Ss 0.996 * 5 0600.13 0.828 ^ 0, 015" 0:034 0.902* 5 027-0058 0.827 ^ 5017 0.046 
rarag [Mpc] 147.287 y 10 3.80 148.277 V oS 3.78 MT.6 55 3 148.8 , 5 $8 


TABLE XVII. 68% and 9596 CL constraints on various free and derived parameters of the PEDE + Ox + M, + Neg scenario 
using several observational datasets. 


stick to this point [19] and the effects of curvature on the 
cosmological models and the key cosmological parame- 
ters should be investigated. Thus, the main theme of 
this work is the curvature of the universe and its impacts 
on extended cosmologies with a special focus on the cos- 


mological tensions. 


other key cosmological parameters. 


In order to investigate these issue, we have consid- 
ered 16 extended cosmological scenarios emerging from 
4 well known cosmological models (ACDM, wCDM, 
wowg CDM, PEDE, see section II for their descrip- 
tions) and performed a systematic study with the use of 
CMB temperature and polarization spectra from Planck 


2018 and its combination with other astronomical data, 
e.g. BAO and Pantheon SNIa catalogue. 


In particu- 


A. Summary 


1. Curvature 


lar, these sixteen extended cosmological scenarios have 
been constrained making use of CMB, CMB+BAO, 
CMB-+Pantheon and CMB4-BAO Pantheon. In the fol- 
lowing we describe the constraints on the curvature and 
then present the effects of the a varying curvature on 


The results for the curvature density parameter are 
summarized in Figure 17 for all the different cosmological 


202302.00238v1 


chinaXiv 


29 


EN PEDE + 9x + Neg: CMB 

Mmm PEDE + O + New: CMB+BAO 

ENH PEDE + Oy + Ne: CMB+Pantheon 

EN PEDE + O5 + Ne: CMB+BAO+Pantheon 


140 - 


L 
2.4 2.8 3.2 3.6 —0.06 


nu Oko Qmo 


0.00 0.30 0.45 0.60 


L l 
0.8 0.9 10 11 48 56 64 72 140 148 156 


Ss H, 0 Tdrag 


FIG. 15. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by PEDE + Qx + Neg using various observational datasets. 


models and data-combinations. 


The Planck indication for a closed universe observed 
within the baseline ACDM case persists in all the non- 
flat extended cosmologies analyzed in this work, although 
at different statistical level because of the large error-bars 
obtained by extending the parameter-space and/or con- 
sidering different parameterizations for the dark energy 
sector. Notice that this loss of constraining power is typi- 
cally due to the strong geometrical degeneracy among dif- 
ferent parameters but in particular among the DE equa- 
tion of state and the curvature density parameter Ngo. 
Such degeneracy can be broken by combining the CMB 
data with other complementary datasets, such as the 
large scale structure information from BAO or the dis- 
tance moduli measurements from the Pantheon sample. 
In the first case, for the combination of CMB and BAO 
(i.e. CMB+BAO), the preference for a closed Universe 
disappears and, regardless from the dark energy model, 


the constraints on curvature shrink around C) = 0, al- 
ways pointing towards spatial flatness (within one stan- 
dard deviation or slightly more for a dynamical dark en- 
ergy). However, it is worth recalling that BAOs are in 
general tension with Planck when curvature is allowed 
to vary and so this data-combination is not robust. On 
the other hand, combining CMB4-Pantheon, the results 
on Qo become sensitive to the assumptions in the dark 
energy sector. Indeed fixing the dark energy equation of 
state to w — —1, this dataset still indicates a spatially flat 
Universe at one standard deviation, while in all the other 
dynamical and non-dynamical dark energy parameteriza- 
tions it systematically shows a preference for a spatially 
closed geometry. 
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FIG. 16. This figure corresponds to the 1-dimensional posterior distributions and the 2-dimensional joint contours for the most 
important parameters of cosmological scenario labeled by PEDE + Ox + M, + Neg using various observational datasets. 


2. Ho 


Figure 18 summarizes the constraints on the Hubble 
constant Ho for all cosmological models and for all data 
combinations showing clearly how different cosmological 
models are returning different values of Hy when the cur- 
vature is taken into account. Throughout the present 
article, we find that the Hubble constant has a deep de- 
generacy with the curvature parameter for the proposed 
cosmological models. 

In all extensions of ACDM, we have an evidence of 
a closed Universe at more than 9596 CL for CMB alone 
dataset. Since Hp and ÎN go are positively correlated in all 
these extended models (see Figure 1, Figure 2, Figure 3 
and Figure 4) and we notice that as long as Qo starts 
deviating from its null value (i.e. Qj = 0) towards the 
negative region, Ho takes lower values due to this corre- 
lation which as a result increases the tension with SHOES 


(see Figure 18). This increased tension is solely due to 
the evidence of a closed Universe for CMB alone in such 
scenarios. For CMB+BAO and CMB+BAO-+ Pantheon, 
the estimated values of Ho in all such scenarios are almost 
same with the Planck value (within flat ACDM), as (ko 
is consistent to zero, with the exception of the PEDE 
parametrization for which Ho is slightly higher alleviat- 
ing the tension with SHOES. For CMB+Pantheon, the Ho 
values in all four extended scenarios are relatively lower 
than the Planck value obtained within the flat ACDM 
scenario. Thus, within these four extended scenarios of 
ACDM, we do not observe any effective resolution of the 
Ho tension. 


Focusing on all the extended scenarios of the wCDM 
cosmology, we find that CMB alone again returns lower 
values of Ho compared to Planck’s estimation (within the 
flat ACDM model) due to an evidence of the closed Uni- 
verse at more than 68% CL (for wCDM + Ox + M,, 
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FIG. 17. Whisker plot with 68% CL constraints on the curvature parameter (Qxo) obtained in various cosmological scenarios 
considered in this article has been displayed for several observational datasets, namely, CMB, CMB+BAO, CMB--Pantheon 
and CMB--BAO-rPantheon. The magenta vertical line corresponds to the flat case (Qo = 0). 


closed Universe is evident at more than 9596 CL). How- 
ever, the error bars in Ho are found to be very large and 
due to this large error bars the tension with the SHOES 
collaboration is reduced down to 1.20 (for wCDM + Ng), 
< 20 (wCDM + Oy + Neg), 2.70 (wCDM + Oy + M, 
+ Neg). However, such alleviation of the Ho tension 
is mainly driven by the volume effect of the parameter 


space. The constraints on Ho from CMB+BAO are al- 
most similar to the Planck values (within flat ACDM), 
however, the mean values of Ho in all the cases and the 
error bars are slightly large than the Planck values and as 
a result we see that for the CMB+BAO analysis, the ten- 
sion with the SHOES for such extended scenarios are re- 
duced down to 2.30 (for wCDM + Ox), 2.40 (for wCDM 


+ Ox + M,), 2.20 (for wCDM + Or + Neg), 2.30 (for 
wCDM + Ox + M, + Neg). For CMB+Pantheon, we 
obtain relatively lower values of Hp than Planck (within 
flat ACDM) due to the preference of a closed Universe 
at more than 9596 CL. Finally, the constraints on Ho 
from CMB+BAO+Pantheon are almost similar to the 
CMB+BAO analysis for all the extended cosmological 
scenarios in this category. 

For wow; CDM + Ox, and its other extensions with 
neutrinos, we see that CMB alone returns lower values of 
Ho compared to Planck's estimation (within flat ACDM) 
but with very large error bars which as a result help to 
reduce the tension with SHOES down to lo, 2.50, « lo, 
2.30, respectively, for wow; CDM + Qg, wow, CDM + 
Ox + ML, wow; CDM + Qg + Neg and wow; CDM + 
Og + My + Neg. For CMB+BAO, in all extended sce- 
narios we obtain a preference for a closed Universe at 
more than 6896 CL and as a consequence we get lower 
values of Ho compared to Planck (within flat ACDM) as 
an effect of this closed geometry. For CMB4-Pantheon, 
the estimated values of Ho are lower than Planck's value 
and they are also lower than the CMB--BAO estimations 
in these extended scenarios. This is due to the fact that 
CMB4-Pantheon supports a closed Universe at more than 
95% CL. For CMB+BAO+Pantheon, the Ho values ob- 
tained in these extended scenarios match almost Planck 
values (within flat ACDM). 

In the extended PEDE scenarios, we see that for both 
CMB alone and CMB+Pantheon, Ho assumes lower val- 
ues similar to the previously explored scenarios and this 
is due to the fact that in both the cases, a closed Uni- 
verse scenario is preferred at more than 9596 CL. For 
CMB+BAO analysis, we see that even though all the 
extended scenarios indicate a preference of a closed Uni- 
verse at more than 6896 CL, however, unlike to the pre- 
vious scenarios where closed Universe preference worsens 
the Hg tension, here on the contrary, the tension with 
SHOES is reduced down to many standard deviations. In 
particular, we find that the Ho tension is reduced down 
to 2c (PEDE + Ox), 20 (PEDE + Ox + M,), 1.90 
(PEDE + Ox + Neg), 1.90 (PEDE + Ox 4- My + Neg). 
The results for CMB--BAO-4-Pantheon are almost similar 
to the CMB+BAO case with a hints of a closed Universe 
model at more than 6896 CL. 

It is worth noticing here that models of new physics at 
low-redshift (or the late time solutions) cannot solve the 
Ho tension once BAO and Pantheon data are considered 
(see Refs. [38, 131, 132]), unless the energy density of DE 
becomes negative [133, 134]. 


3. Ss 


In Figure 19 we have summarized the constraints 
on the Sg parameter extracted out of all the ex- 
tended cosmological models using CMB, CMB+BAO, 
CMB+Pantheon and CMB+BAO+Pantheon together 
with the estimated values of Sg from Kilo-Degree Survey 
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(KIDS-1000) [Sg = 0.766*0 074] [71], Dark Energy Survey 
(DES) Year 3 data (DES-Y3) [Sg = 0.759*0025] [73] and 
Planck Sg = 0.834 + 0.016 [8], all obtained by assuming 
a ACDM scenario as the underlying model. 

From Figure 19 we observe two things. For CMB and 
CMB-+Pantheon, the estimated values of the Sg param- 
eter in all the extended cosmological scenarios are very 
far from both KIDS-1000 [71] and DES-Y3 [73], when 
they assume a vanilla ACDM model.? Moreover, they 
are also far from Planck's estimation assuming a vanilla 
ACDM model, except in the extended scenarios assum- 
ing a cosmological constant, where for CMB+Pantheon 
case Sg values are slightly close to it. For all the other 
cases, for these two datasets the tension in the Sg pa- 
rameter increases significantly. This is because of the 
indication for a closed Universe coming from these two 
datasets when it is not assumed a cosmological constant. 
On the other hand, when we look at the estimations of Sg 
by CMB+BAO and CMB+BAO-+Pantheon, we see that 
these estimations are almost similar to the Planck values 
assuming a vanilla ACDM model. This indicates that 
within such extended scenarios, the Sg tension increases 
when there is an indication for a closed Universe. 


4. M, 


From all the analyses, we find no evidence of a to- 
tal neutrino mass. In particular, only in a few cases, 
namely, ACDM + Qs + M, + Neg, wCDM+ Oy + M,, 
wow; CDM + Ox + M,, PEDE + Or + My + Neg, for 
CMB alone, we find an indication of a total neutrino mass 
at 6896 CL but this evidence disappears when other exter- 
nal and complementary probes are added to CMB. How- 
ever, in order to realize the bounds on M, in all these sce- 
narios, in Figure 20, we have presented a graphical pre- 
sentation of the the upper limits of M, at 9596 CL for our 
extended cosmological models using CMB, CMB+BAO, 
CMB-+Pantheon and CMB+BAO-+Pantheon. 


5. Neg 


We provide a summary picture of Neg in Fig- 
ure 21 obtained in all our extended cosmological sce- 
narios considered in this work for CMB, CMB+BAO, 
CMB-+Pantheon, and CMB+BAO+Pantheon. We find 
that irrespective of all the cosmological models con- 
sidered in this work, the estimated values of Neg for 
CMB and CMB+Pantheon are almost same with that 
of the standard value Neg = 3.044 and thus these 


6 Both KiDS-100 and DES collaborations explored extended sce- 
narios in Refs. [72, 74], but the combinations of parameters ex- 
plored here are not considered, so we can't make a safe compar- 
ison. 
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FIG. 18. Whisker plot with 6896 CL constraints on the Hubble constant Ho for a variety of cosmological scenarios in 
the presence of a non-zero curvature of the universe has been shown considering CMB, CMB+BAO, CMB+Pantheon and 
CMB+BAO+Pantheon. The yellow vertical dotted line corresponds to the Ho value from Planck 2018 team [8] assuming a 
ACDM model and the magenta vertical band corresponds to the Ho value from SHOES team (Ho = 73.044 0.99 km/s/Mpc at 


6896 CL) [32]. 


estimations are consistent with this standard value. 
For the remaining two cases, namely, CMB+BAO and 
CMB+BAO-+Pantheon, the values of Neg are slightly 
lower (Neg < 3) but within 68% CL, they are also con- 


sistent with the standard value. 


6. wo and wa 


Here we summarize the main results on the dark energy 
equation of state parameter w of the wCDM model and 
also on the (wo, wa) parameters of the wowa CDM model. 
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FIG. 19. Whisker plot with 68% CL constraints on the Sg (= 08\/Qmo/0.3) parameter obtained in various cosmological scenarios 
considered in this article has been displayed for several observational datasets, namely, CMB, CMB+BAO, CMB+Pantheon and 
CMB+BAO+Pantheon. The magenta vertical band corresponds to the estimated value by the Kilo-Degree Survey (KIDS-1000) 
leading to Sg = 0.766+9:929 [71], while the green vertical band corresponds to the value provided by the Planck collaboration 
Ss = 0.834 + 0.016 [8], all obtained by assuming a ACDM scenario. 


In Figure 22, we show the whisker plots with 6896 CL, 
and 9596 upper limits, on w and wa for the wCDM and 
voa CDM models for all the observational datasets. 


One can clearly see that in the wCDM model, only the 
CMB-+Pantheon dataset gives a strong indication of a 


phantom DE (i.e. w « —1) at more than 9596 CL, only 
when the equation of state of the DE is constant and not 
dynamical. For the remaining cases, w — —1 is allowed 
within 6896 CL. 


However, for the wow; CDM model in which the dark 
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FIG. 20. Whisker plot with 95% CL upper bounds on M, [eV] obtained in various cosmological scenarios considered 
in this article has been displayed for several observational datasets, namely, CMB, CMB--BAO, CMB-Pantheon and 
CMB+BAO-+Pantheon. 


energy equation of state is dynamical, CMB+BAO gives ^ with a cosmological constant’. Regarding the dynamical 
an indication for a quintessence DE (ie. wo > —1), nature of the DE quantified through the free parame- 
while all the other dataset combinations are in agreement 


7 Note here that wo remains unconstrained in the scenarios 


€ CMB 
B] CMB-BAO 


@ CMB+Pantheon 
@ CMB+BAO+Pantheon 


36 


. 


PEDE--Oy + M, + Nest 


3 d 


PEDE+Qx + Nes 


i 


WoWaCDM + Qk + My + Neff 


WoWaCDM + Ox + Ner 


E —————4————————À4 
p — ee 

co 
IM 8—— — — —3À4 

e 
—— e —À 


41 ! 
I 


—— —4.8— —— ————3À 


wCDM + Qk + M, + Ner 


wCDM + Ox + Ner 


ACDM + Oy + My + Ner 


-j-X-r---9- 


l 


ACDM + Ox + Neff 


| 


FIG. 21. Whisker plot with 68% CL constraints on Neg obtained in various cosmological scenarios considered in this article has 
been displayed for several observational datasets, namely, CMB, CMB+BAO, CMB+Pantheon and CMB+BAO-+Pantheon. 
The magenta vertical line corresponds to the standard value Neg = 3.044. 


ter Wa, we find that CMB alone and CMB+Pantheon 
are very consistent with wa = 0 (except in the model 
woWaCDM + Qx in which wa z 0 at more than 68% 
CL), however, when BAO data are included to CMB, 
an indication for a wa < 0 at more than 2c is found. 


woWaCDM + Or + My, wowa CDM + OK + My + Neg and 
for the scenario woWaCDM + Ng + Neg, even though the 68% 
CL upper bound on wo is available but it becomes unconstrained 
at 95% CL. 


While this evidence is mildly diluted (but remains true 
at slightly more than 68% CL) when the full combina- 
tion CMB+BAO-+Pantheon is explored. In other words, 
for this wow; CDM model, CMB+BAO is preferring a 
quintessence dynamical DE at more than 95% CL. 


B. Concluding remarks and the road ahead 


There is no doubt that without observational data, 
modern cosmology is incomplete. From the detection of 
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FIG. 22. Whisker plot with 68% CL constraints and 95% CL limits on wo and wa obtained in various cosmological scenarios 
considered in this article has been displayed for several observational datasets, namely, CMB, CMB+BAO, CMB+Pantheon 
and CMB+BAO+Pantheon. The magenta vertical lines in the two panels represent to the standard cosmological constant 
value wo = —1 and to the non dynamical wa = 0 limit, respectively. The missing blue lines (Planck only results) in the left 


panel are indicating that w is unconstrained for those models. 


the cosmic microwave background radiation to the dis- 
covery of the late-time cosmic acceleration, we have wit- 
nessed how the observational data have been crucially 
important in the understanding of our Universe. When- 
ever we talk about cosmology, the six parameter ACDM 
cosmological model under the assumptions of a positive 
cosmological constant, a pressure-less DM and a flat Uni- 
verse, naturally enters into the picture because of its su- 
per agreement with a series of astronomical data. How- 
ever, the true nature of DE and DM within this simplest 
cosmological scenario has remained one of the longstand- 
ing debates in cosmology. In addition to that, the as- 
sumption of a flat Universe has created a further debate 
in the scientific community. A series of recent articles 
have argued that the present observational data hint to- 
wards an evidence of a closed Universe [8-11, 18]. This 
evidence has sparked the entire scientific community be- 
cause any evidence of a closed Universe increases the ten- 
sions in the key cosmological parameters and hence the 
flat ACDM cosmology is further challenged. Generally, 
the curvature of the Universe is “assumed” to be zero 
when the cosmological models are analyzed. This is mo- 
tivated because: (i) the inflationary paradigm has been 
tremendously successful in demanding a spatially flat 
Universe, (ii) the observational data in the past preferred 
a small curvature of the Universe which is very close to 


zero, and further (iii) the inclusion of extra parameters, 
in terms of the curvature parameter for instance, may in- 
crease the degeneracies in the parameters. However, over 
the last several years, the sensitivity of the astronomical 
data has improved a lot, and given such improvements, 
it is very natural to allow the observational data to de- 
cide the nature of the curvature parameter. That means, 
the consideration of a flat Universe (i.e. Qg = 0) while 
analyzing the cosmological models is not a very realistic 
assumption. We further comment that forcing Qg = 0 
could bias the results and the intrinsic nature of the cos- 
mological models may remain hidden forever [19], and at 
this stage where we have a large amount of astronom- 
ical data from various sources, having Qg = 0 is not 
expected. 


Thus, in the present article we investigated the effects 
of curvature on the cosmological models and their associ- 
ated (free and derived) parameters. We considered some 
classical and well known cosmological models and their 
several extensions, including variations in the DE and 
neutrino sectors, in presence of the curvature of our Uni- 
verse. Our analyses clearly indicate that it is better not 
to fix Or = 0, to analyze the cosmological models. We 
trust that the upcoming observational data from various 
astronomical surveys (see for instance, Refs. [135-138]) 
will be highly promising in this direction which will re- 
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veal the intrinsic nature of the dark sector as well as the 
curvature of the Universe. 
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